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I . INTRODUCTION 


Satellite laser ranging systems are currently being used to accurately 
measure the baseline distances between widely apearated points on the 
earth's surface [I]« The technique involves measuring the distances be* 
tween the ground stations and retroref lector equipped satellites [2]. 
Typically, hundreds of range measurements are obtained during a single 
satellite pass over a ground station* Measurements from several passes are 
used to determine the satellite orbit and the coordinates of the ground 
sites* Major error sources include Instrumentation noise, orbit modeling 
errors and atmospheric refraction* The effects of these errors on baseline 
accuracy is a complicated function of the geometry of the satellite orbits 
and ground station locations* In this paper we are primarily concerned 
with the effects of atmospheric refraction. Unfortunately, because of the 
mathematical complexities involved in exact analyses of baseline errors, it 
is not easy to isolate the atmospheric refraction effects. However, by 
making certain simplifying assuL^ptions about the ranging system geometry, 
relatively simple expressions can be derived which relate the baseline 
errors directly to the refraction errors. The results indicate that even 
in the absence of other errors, the baseline error for intercontinental 
baselines can be more than an order of magnitude larger than the refraction 


error* 


2 


2. BASELINE ERROR ANALYSIS 

The general nonlinear regression model of n measurements 13] is given 
by 

i ■ ^ ± ^ ^ ^ 

where jc contains epoch values of the desired parameters to be estimated, 
contains epoch values of the unadjusted parameters which are assumed to be 
known constants in solving the regression equations, and ^ is the zero mean 
measurement noise vector. It is assumed that the elements of ^ are sta* 
tistically independent, and the partitioning of parameters into and ^ is 
arbitrary. 

When this regression equation is used to model the distance measure* 
meats between a ground station and the satellite, X is a vector repre- 
senting the range measurements between a ground station and the satellite, 
X is a vector representing the epoch values of the ground station coor- 
dinates, and is a vector representing the epoch values of the error 
sources' parameters such as the orbit modeling errors and the atmospheric 
refraction errors. The vector function ^( 3 ^,^) in Equation (1) can be 
defined as 

where ^ is a vector representing the r.ctual geometric distances between the 
ground stations and the satellite, and ^ is a vector representing the 
atmospheric refraction. 

The total error covariance matrix associated with the station coor- 
dinates for this model is given by 131 
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EKAxXAx)’’^] - (a’W+ + (A^WA+ (bJwA) 

T -1 -1 

«(A WA + P ) 


(3) 


where 


3f(x,s) 



s 




(4) 


B 


3s 






(5) 


Is the covariance matrix for the a priori estimate of the station coor- 
dinates ^ is the inverse of the measurement noise covariance matrix, 

V is the covariance matrix associated with the unadjusted parameter s, and 
— s — 

and ara the nominal values for x and s, respectively* 

— ^ — N — — 

The elements of the matrix V[ are directly proportional to the total 
number of range measurements. So in satellite laser ranging where a large 
number of range measurements is obtained, the matarix in Equation (3) 
can be neglected when it is compared to the matrix. That is, the a 

priori estimate of is not as critical in determining the error covariance 
matrix when a large number of range measurements is gathered. Because 
of this and the fact that the (A WA) matrix is smal'. , Equation (3) can be 


simplified to 
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( 6 ) 


In Cartnuian coordinates, the baseline distance between two points can 
be expressed as 

\ ’ ^ 2 ^^ “*■ ^*1 ” 

where (x^,y^,*,) and coordinates of the two points of 

interest. By letting 4x , 4y , Az denote the coordinate measurement 

tft fli (II 

errors associated with the mth station (m ■ 1,2), the baseline error can be 
expressed mathematically by expanding Equation (7) in a Taylor series about 
the current nominal parameter x^. 


B » — ^ 
L 3Xj 




* -is 

i’is 


Ayi 


i- *N 





+ Si 

^2 







AZ. 


( 8 ) 


i'iSN 


It can be seen from Equation (8) that the mean square baseline error 
consists of the variances and the covariances of the various coordinate 
errors. By explicitly evaluating Equation (6) and using the results with 
Equation (8), the baseline error can be determined. 
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3. RANGING GEOMETRY 

As asntlonsd in Section 1, the osjor error sources for the bssellne 
deteralnetlon Include the orbit oodellng errors end the staospheric refrac- 
tion errors. In this analysis, we want to reduce all other error sources 
to negligible levels in order to estimate the ultimate effect of refraction 
on the baseline errors. By doing so, we are assuming that the orbit modeling 
errors are negligible. This Is a reasonable assumption because. In practice, 
several days' worth of data from many stations are used to determine the best 
fit satellite orbit (4). The station coordinates are then determined using 
this best fit orbit. In this section, a particular ranging geometry for the 
best fit satellite orbit Is defined and used to calculate the baseline error. 

The geometry of the two ground stations and the two satellite paths Is 
illustrated In Figure 1. Notice that two satellite passes are required to 
solve for the ground station coordinates. The satellite ground tracks are 
parallel to each other, and the satellite altitudes (h) are constant and 
equal for the two passes. This corresponds roughly to the geometry for a 
polar, circular orbit satellite. Furthermore, high altitude orbits are 
assumed so that a flat-earth model can be employed In the analysis 
to simplify the mathematics. It should be noted that the results otalned 
by using the flat-earth model will not differ substantially from those 
obtained by using the spherical-earth model. 
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Figure 1. Geometry of the two-ground-statlon, two-satelLlte-pass ranging 
system. 
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For slaplicicy» ch« aactlllt* orbict, and «r« Msuatd to b« in 
cho x*dlr«cclon. Th« noaliul coordlnatao of ch« ground sctclon 1 (C^) ar« 
((L^ * L^} ten Yi 0) and the noalnal coordlnacat of th« ground acaelon 
2 (G^) are (O.'L^iO)* The baaallna batwaan tha cao acatlona and cha y*axia 
Incarsacc at an angle y* ^a aaxlaum alavaclon angle asaoclatad with Cha 
fflCh ground staclon and cha nch aacalllca pass la danocad by 

a - 1.2. n - 1,2 . 

nn 

For chla parcicular gaoaacry , 

ix ■ (-1)**^^ dx ain Y a ■ 1.2 , (9) 

m a 

Ay ■ (-1)*^^ Ay coa y a ■ 1,2 , (10) 

n ® 

and 




AT - 0 a • 1,2 . (11) 

a 

i-i£N 

By subscicuclng EquaClona (9) through (11) inco Equation (8), we have 
AB^ » sin Y (AXj - Ax^) + coa y (Ay^ - Ay^) . (12) 



The mean square baseline error la, therefore, given by 
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4. BASELINE ERROR FOR THE SPHERICALLY SYMMETRIC ATMOSPHERE 
The various matrices m^ntloaed In Section 2 will be formulated here 
for the ranging geometry of Section 3. For this ranging model, the orbit 
modeling errors have been assumed to be negligible and the dominant errors 
associated with the range measurements are the atmospheric refraction errors 
and the instrumentation errors. 

Normally, the range measurements are taken when the elevation angles 

are above certain threshold angles E^ and they are taken uniformly in time 

inn 

during the satellite pass (refer to Figure 2). Thus, we can denote the 
position of the satellite during the jth measurement by its Cartesian coor- 
dinates as 

if 1 < j < k 

^2 

if k+1 < J < 2k 
S 

I 

if 2k+l < J < 3k 

^2 

if 3k+l < j < 4k. 

It should be noted that the threshold elevation angles associated with each 
ground station and each satellite pass are different, but they are usually 
related to one another for a particular geometry. Further discussion on 
their relationship is given later In this section. 


position of satellite on orbit 
during the jth measurement for 
station 




position of satellite on orbit 
during the Jth measurement for 
station 

position of satellite on orbit 
during the jth measurement for 
station G. 


position of satellite on orbit 
during the jth measurement for 
station G. 
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If Gj and are located at (x^.y^z^) and (X2,y2**2^» r«*P«ctively, 
then the geometric distances between the ground stations and the satellite 
at various points along the satellite orbits can be written as 


[(x_.-x,)^-»'(y„.-'y,)^+(h-z,)^l^ geometric distances between G. and 
^ ^ ^ satellite along if 1 < J < K 

l(x -x,)^+(y-.-y,)^+(h-z,)^l^ geometric distance between G. and 
^ ^ satellite along if k+l < J < 2k 


l(x_.-x,)^+(y„.-y,)^+(h-z»)^l^ geometric distance between G, and 
^ ^ satellite along if 2k+l < f < 3k 


l(x -x_)^+(y-.-y,)^+(h-z-)^)^ geometric distance between G and 
^ ^ ^ satellite along if 3k+l < J < 4k, 

(14) 


and the column vector ^ of Equation (2) can be expressed as 

I d|^ , . . , • (15) 

For a spherically symmetric atmosphere, the effects of horizontal 
refractivity gradients 1 3] are negligible, and Che atmospheric refraction 
associated with Che jch range measurement is given by 16] 

AR i — j ■ l,...,4k (16) 

^ sin Ej 

where E^ is the elevation angle associated with the ground station and the 
satellite at the jth measurement, and Bj fs the spherical correction coef- 
ficient associated with the jch measurement. Bj is a function of surface 
pressure, temperature and Che water vapor pressure at the ranging sites 
during Che satellite passes [6]. By referring to the ranging geometry of 
Figure 2, Equation (16) can also be expressed as 
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j ■ 1 , . . . ,4k . 


( 17 ) 


By using Che results of Equation (17), the column vector^ of 
Equation (2) can now be expressed explicitly as 


AR - 


0.d. 




ri 


’4k"4k 

TJ 


(18) 


By substituting Equations (15) and (18) into Equation (2), the range 
distance vector spherically synaetric atmosphere can be writ- 

ten as 


f (x,s) 



(19) 


The spherical correction coefficients vary from measurement to 
measurement, depending on the instantaneous weather conditions at Che 
ranging sites and along Che optical path. So Che unadjusted parameter vec- 
tor 3 should be expressed as 





( 20 ) 


Since the tine period for each satellite pass Is usually shorty substantial 
changes in the weather ^.onditlons at the ground stations are not expected. 
Typically » meteorological data Is acquired once during each satellite pass 
and used to calculate B to correct the ri*aging data collected during that 
pass. Therefore 9 ve assume that the value of B that is used to correct 
the measurements for a ground station remains constant on the same 
satellite pass. 
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The ground stations are usually widely separated so that the surface 
pressure, temperature and the water vapor pressure at these stations eay be 
coepletely different froe and Independent of each other. As a result, the 
aeteorologlcal data, and hence the 0's, associated %rith different ground 
stations, are assuaed to be Independent. Consequently, the errors In calcu- 
lating these 3's (A 0's) at the two ground stations are assuaed to be 
uncorrelated. Notice that these errors arise froa Inaccuracies In the for- 
auias used to predict 0 froa aeteorologlcal paraiseters and froa errors In 
the meteorological aeasurenents. 

Due Co Che assuaptlons of constant 0 for Che saae satelilCe pass and 
dlfferenc 0's becween ground scatlons, Equaclon (20) can be simplified as 

s.- 10^^ 0^2 821 

where Q is the estimated spherical correction coefficient associated with 
mn 

the mth ground station during the nth pass. 

Because the range measurements in the two passes are taken at dif- 
ferent times, the correlation of A0 between passes can either be high or 
low, depending on the extent of the weather changes during the two passes. 
When the temporal separation between passes is short or when the weather 
conditions at the ground stations are fairly constant during the two 
passes, we would expect to see a very high correlation betwen the values of 
B for the two passes. Therefore, the errors in B» i.e*, AB, would be 
expected to be correlated for the two passes. Conversely, if the time 
separating the two passes is long or if weather conditions changed substan- 
tially, a partial or a relatively low correlation between AB during the two 
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passes would be expected. In Section 4.1, a general expression for the 
baseline error In the spherically syiaeecrlc ataosphere Is derived. 

4.1 General Expression for the Baseline Error 

Before we can derive a general expression for the baseline error, the 
various matrices of Equation (3) must be expressed explicitly in terms of 
the relative positions of the satellite and the ground stations according 
to Che prescribed geometry. 

By replacing the general B's of Equation (20) with Che station-varying 
and path-varying 6 's, we obtain a simplified expression for Che range 

IBll 

distance vector as 


f(x,s) 





J 



1 + 



r 

1 + 



1 + — ,... , 

®2ll 

1 + 

. 

d- . 

f ®22l 


®22l 

1 + 

d. 

h ! 2k+r 

h , 

L / 

3k 

h 

®3k+l ’ * * * ’ 

^ J 

**4k. 


( 22 ) 


where d^ is a function of both Che station coordinates and the current 
satellite positions, and 

X- Ui ' 2 '^ (23) 

is the ground station coordinate vector. 

The partial derivative of with respect to jc can be calculated 

by differentiating each vector component of Equation (22) with respect to 
each vector component of Equation (23). For instance , differentiating the 
first component of Equation (22) with respect to the x-coordinate of the 
ground station 1 yields 
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31U +-^) djl 

3x, 


1 + 


’ll 


h J 


(24) 


3f 


Notice that the resulting partial derivative is a 4k by 6 matrix. 

After evaluating this matrix at the nominal station coordinates of 
((L^ L2> tan y. 0) and (0,-L2,0)(refer to Figure 1), the ^ matrix 

of Equation (4) can be expressed explicitly as 



(25) 


The partial derivative of with respect .o ^ can be calculated by 

di .ferentiating each vector component of Equation (22) with respect to each 
vector component of Equation (21). The resulting matrix is a 4k by 4 
matrix. After evaluating this matrix at the nominal station coordinates, 
Che B matrix of Equation (5) can be expressed explicit Ly as 





0 ... 0 


\-K ^ 

h *•• h 


21cfl 


*3k 


3k+l 


h 


16 


(26) 


In order to derive an explicit expression for the natrix, we have to 
Investigate the statistics of the measurement noise in the satellite laser 
ranging system. It has been shown that the phctoelectron distribution from 
a single-mode laser follows Poisson [7]. Recently, the MMSE and ML estimator 
noise variances for the laser ranging receiver have also been derived [8]. 

In particular, if the received optical pulse shape closely approximates a 
Gaussian pulse shape, and if no background noise is present in the system, 
the ML estimator noise variance for each range measurement is given by [8] 



(27) 


where b is the full width of the optical pulse measured at half maximum 
(FWHM), and Q is the average number of photoelectrons per pulse. 

The average number of photoelectrons in a received pulse is directly 
proportional to the received pulse strength, which in turn Is Inversely 
proportional to the geometric distance between the ground station and the 
satellite. In a vacuum, this signal pulse strength decreases inversely 




i 


- — - 
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wlch ch« squara of Cha distaaca. In tha aarth's acaosphara whara abaorp* 
tion and scattaring incraasa with pathlangth, tha signal pulsa strangth nay 
dacraase avan mora quickly. 

In tha satalllta lasar ranging systaa, both the laser transalttar and 
the receiver are located at the ground stations, and ratroraf lectors are 
Installed on tha satalllta. The transalttar sands out pulses at a constant 
repetition rate during the ranging period. These pulses are than reflected 
back from the retroref lector-equipped satellite and are detected by the 
receiver. The tlae of flight of the laser pulse Is measured through Che 
range counter, which Is Chen mulclplled by Che velocity of light in the 
aCaosphere Co give the range to Che satellite [2], [9]. 

In this configuration, Che optical pulses travel twice Che sCatlon-to- 
satelllte dlsCat'ce. So Che received pulse strength, and therefore the 
number of photoelectrons per pulse. Is Inversely proportional to the fourth 
power of Che distance, and the measurement noise variance can be expressed 
as 


Oj - dj J - 1 4k (28) 

2 

where o^ Is called Che aaasureaenc noise variance factor. By using the 
result of Equation (28), we can obtain the Inverse of the measurement noise 
covariance matrix as 


dlag 

a. 




(29) 


By the previous assuapClons that Che errors in B's (46) between sta- 
tions are uncorrelated and that the degree of S*error correlation between 
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passM d«p«adt upon diffcrcncas la Che sceciooe' weather condlcloae during 
Che two paesee , Che covariance nacrix assoclaced with Che vector ^ can be 
expressed as 


‘ 0 
0 : 


(30) 


where ot is Che variance of 8, and p , m ■ 1,2 is Che 8-error correlacion 

P 1ft 

coefficient associated with Che aeasureaents acquired for Che mch station 
during Che two passes. 

By subsClCuCing Equations (25), (26), (29) and (30) into Equation (6) 
and carrying out the appropriate satrix operations, the error covariance 
aacrix associated with the a cation coordinates can be obtained as 


El(^)(^)^J • 


cov(yj,*^) 


cov(y^,z^) 

0 

0 


0 

0 


0 

0 

0 

0 

0 


0 

0 

0 

0 


cov(y2,*2) 


0 

0 

0 

0 

cov(y2,Z2) 

2 

a 


(31) 


where 


.1.2 e'Jj .1.2 

, 2,_M pM V 

sin (E^^ - Ej 2 > 


’ll 


2piCiiCi2 


4 M 

Sin Ej^ 


sin^ e”j sin^ E^^ 


'12 


sin'* E ‘^2 


,( 32 ) 
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-•ovCy^.z^) 


o; .in‘ E^^ »la‘ E”2 


PjCjiCj2 

- ‘u> 

•in^ eJj sin^ E^2 

tin^ eJ^ sln^ E^2 


•(cot + cot E^,) + 


'12 


sin^ E‘^2 


(33) 


2 , 2 -M , 2 pM 

Og sin E^ sin E^2 


'll 




sln^ eJj tan^ 






sln^ E^^ sin^ E^2 ®i2 




sin^ E‘^2 


(34) 


2 , 2 _M 

agsin E^^ 

.1,2 eJj 


2O2C21C22 

^ ^22 

sln2(E”^ - £^2^ 

.ln‘ E« 

sln^ E^j sin^ 

E^2 .ln^E^2 


,( 35 ) 
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cov(y2,*2) “ 


4 b;, 

, 2,,M M , 
•In (E 2 i“ ®22^ 


*21 


•in ^21 ^22 


PaOa • 0^ M 
2 21 22 

. 2 _M , i _M 
•In E 2 ^ >ln E 22 


•(cot + cot £^ 2 ) 


'22 


_M , ^ 

•in E 22 tan £ 2 ^ 


(36) 


2 



-2 _M ..2 „M 

^ ^21 ^ ^22 

•ln^(E 2 i - E22) 


— , 
0 


2: 


4 2 

M 

ain E 2 ^tan 

®22 


2 1 ^<1 A 

2 21 22 


, 2 „M ^ 2 ,M ,M 
•in E^^ tin E^^ tan E^^ tan E^^ 


'JL 


, 4 „M ^ 2 _M 

sin B 22 tan 


(37) 


c T‘^° *2 ».n 

mn 1 j j 


• •1.2, n-1.2 (■;«) 


and 


’^an 


oo« 


-1 


sin E 


an 


sin E 


M 

'an 


a ■ i ,2, n • 1 ,2 


(39) 


Equation (38) !• plotted versus ^ in Figure 3. This figure shows 

nn 

that Che effective hounds for C are 

an 

I < ^ ‘ •T radians . 

sn tun 2 
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Flgurt 3. The effective bounde for C^. 
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Thl« sMll r«ng« of C indicotoo that C ic vory to ehongoo 

lUi n 


of vorlotto poraaottro. 

Equation (31) indieataa that tha poaltion arrora batwaan tha two 

ground atationa ara uncorralatad. Thla la to ba expactad alnea tha 

■ataorologlcal data arrora and aaaaurwMnt arrora ara uncorralatad batwaan 

tha wldaly aaparatad atationa* Furthamora, baeauaa of sysAatry, tha 

x-coordinata arrora nada on tha two paaaaa off anting aach othar, and tha nat 

x-coordinata ranga arror (tha x^axia ia parallal to tha aatallita patha) 

2 2 

for aach atation ia aaro* Tharafora, , cov(xpy^)t cov(xpt^}, o^, 

cov(x2,y2) And eov(x2,t2) all hava xaro antriaa. 

By aubatituting Bquationa (31), (32), and (35) into Equation (13), an 

2 

axplicit axpraaaion for <dB^> for tha aatallita laaar ranging ayatan in tha 
apharically aynaMtric atnoaphara can ba obtainad 




2 2 
Og cos r 


ain^ E?? 


* 2 _M 
11 ®12 


aln^(Eyj - Ej2> 


'11 


2piCiiCi2 


ain^ Ejj 


ain^ E^^ ain^ E^2 


''12 


a in E 


12 


, 2 -M , 2 -M 
ain E2^ ain E22 

•in (Ej, - Ejj) 


'21 


ain 




21 


2P2C21C22 

. i _M • 2 _M 
ain E2 j ain E22 
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4 4 _M 
•in Ejj 


(* 0 ) 


In ganaral, tha aaxlnua alavation angla ia not tha apacifiad paraaatar 


in tha aetalllca laaar ranging procaaa. Tharafora, Equation (40) auat ba 
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sodifiad CO becoM a function of «om known pnrnatcnrs luch n« the minim* 
elevation angle and Che orbit aeparaClon before we can exaaine its proper- 
ties* To illustrate this, Equation (40) is evaluated for the special case 
where the two ground stations are equidistant froa the x-axis* 

In this special case, let the satellite orbits be i/2 away froe Che 
y • 0 plane, and the two ground stations be equidistant froa the x-axls. 
The baseline distance is This geoaetry is illustrated in Figure 4. 

In order Co acquire the saae nuaber of range neasureaents on the two 
passes, Che threshold elevation angles between these passes aust be dif- 
ferent, i.e., * E^. By referring to Figure 2, the specified ainlaua 

elevation angle of Che ranging geoaetry aust be equal to the 

threshold elevation angle formed by a station and a satellite path farther 
away froa that station. That is. 




(41) 


The reaaining threshold elevation angles can be found as 


E°2 ■ - tan"^ 


h Can E 


Bin 


(h^ - CO. , 


(42) 


The sine of Che various maximua elevation angles can be expressed as 


4 4 

sin - sin E 22 


[h^ + V4(Bj^ cos Y - i)^J^ 


(43) 


and 


sin E ^2 - ‘TJ 


(h^ + V4(Bj^ cos Y - i)^j^2 


(44) 


2t 
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Figure 4. Ranging gaoaetry for the case where the two ground stations are 
equidistant from the x-axls. 
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By using ths results of Equations (41) through (44), the equalities 


sin^ sin^ E^2 ®2l ^22 


8in^(E5'j - Ej2) 


sin^(Ejj - E^2> 


( 45 ) 


and 


( 46 ) 


Cii - C22 , Cj 2 - C21 . 

caa be derived# 

By subecituting Equations (43) through (46) into Equation (40) and 
taking the square roots on both sides, an expression for the nns baseline 
error can be obtained as 

2 l/i *^2 Ofl cos Y. 


\ - < 4 Bl> 


'0 

"nr 


(CjJlh^ + V4(Bj^ cos Y + 1)^1 


2,2 


- (pj + P 2 ) Cj^C^ 2 lh^ + l/4(Bj^ cos Y + D^llh^ + I/4(Bj^ cos y “ 0^1 

In order to relate the baseline error directly to the atmospheric 
refraction error. Equation (47) is normalized by the maximum atmospheric 


refraction error (o^ 
elevation angle). This gives 


Bin 


B, /2 cos Y sio , 

L _ Bin j-2 ,.2 


, the refraction error for the ainiBua 


2,2 


AR 


ih 


{Cjjlh^ + 1/4(B, cos Y ”*• « J 


- (pj * P 2 ) C^jCj 2 lh^ * V4(Bj^ cos Y + D^Jlh^ + l/4(Bj^ cos y * 0^1 
+ Cj 2 lh^ + l/4(Bj^ cos Y - i)^l^}^ 
where is a diaensionless parameter called the baseline error 

'’ar 


26 


■ttlciplltr. Wh«n th« Mxiaua r«fr«ctioo «rror aultlpll«d by this 

■ultiplltr, th« aagnlcuds of tho rma botollno orror con b« obtained • 

Equation (68) indicataa that tha baaalina error ia aaallar for the 
higher correlation of refraction arrora on tha two paaaaa, and vice varaa* 
In all caaaa. tha baaalina arror ia an Incraaaing function of tha baaalina 
dtatanca, and la proportional to coa y. 

Tha Incraaaing charactarlatic of tha baaalina arror with raapact to tha 
baaalina diatanca can ba axplainad by obaarvlng that tha Incraaaa In baaa- 
lina distance cauaas an Increase In tha statlon-to-satallita distancas. 
Increasing tha statlon-to-aatalllta distancas will result In larger range 
aeasureaent errors. Consequently, tha baseline error, trhich depends solely 
on the range aeasureaent errors, will also be Increased. 

The dependence of baseline error on angle y can be explained by 
observing that the baseline error Includes only tha y-coordlnata errors. 

Tha 2-coordlnate errors contribute second-order effects which are not 
Included In this analysis. By keeping the other paraaatars constant, tha 
closer y approaches 0 *, the longer tha y-coordinata baseline will ba, 
which will result In a larger baseline error. When y equals 90*, the y- 
coordlnata baseline will be sero. As a result, the corresponding base- 
line error will be sero. 

An approxlaate lower bound for the baseline error can be found by 
assualng a zero baseline and letting be unity. This yields 


> 1 ^ cos Y sin - (o^ + P 2 )l 

®AR 


1/2 I" h 



(49) 


Inequality (49) sho%ra tha dapandanca of baaalina arror on tha orbit 
gaoaatry axplicitly. It also indicataa that a baaalina arror of aero for 
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an arbitrary gaoaatry can only ba achlavad %rhan Cha refraction arrort on tha 
two paaaaa ara parfactly corralatad. 

4.2 Baaallna Error for the Uncorralatad-Path Laaar Ranging 

Whan tha ground stations have undergone large fluctuations In surface 
pressure, twperature and water vapor pressure on the two satellite passes, 
a sero correlation of aeteorologlcal data errors on the two passes will be 
expected. Consequently, the ataospheric refraction error on the two passes 
will also be uncorrelated. This Inplies that 


p • 0 a • 1 ,2 . (50) 

IB 

We refer to this configuration as the uncorrelated-path ranging. The nor- 
aallzed ras baseline error for this configuration can be obtained by 
substituting Equation (50) Into Equation (48). This yields 


% 


(uncorr) 


AR 


7T cos Y slQ 

5i 


— (cj^lh^ + l/4(Bj^ cos Y + 1)^1^ 


Cj 2 lh^ + V4(Bj^ cos Y - 1)^1^}^^ • 


(51) 


By substituting Equation (30) into Inequality (49), we can express the 
lower bound for the baseline error in the uncorrelated-path ranging system 
as 


— ^ (uncorr) 
®AR 


> 2 cos Y sin E . 

aln 



(52) 


Inequality (52) indicates that a nonzero baseline error Is always pre- 
sent (except when y equals 90* ) no aatter how saall the baseline is. 

The baseline distance is liaited by the fact that both satellite passes 
Bust be visible at both ground stations. Consequently, the baseline 
distance cannot be exceeded by Che value 
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B 


Lmx 


2h_ 

tan E 


- t 


aln 


cos Y 


(53) 


With the result of Equation (53), we can find the upper bound for the base- 
line error as 


-B, 


i (uncocr) < I . 


(54) 


31 sin E 


"AR min 

The nomalized ras baseline errors for the uncorrelated-path laser 

ranging system versus the baseline distance are plotted in Figures 5 

through 13 for different sets of paraoeters. As expected, these plots 

Indicate that the baseline error is an increasing function of the baseline 

distance; and it is a cosine function of the angle y* 

The exact dependence of baseline error on the orbit separation and the 

orbit altitude is very difficult to see in general. But for the short 

baseline ranging, the baseline error is shown by Inequality (32) to be 


il + 

1 4h 


, and the oiniaua base- 


related to the two paraneters according to 
line error can be achieved when the orbit separation is twice as ouch as 
Che orbit altitude, i.e., 1 • 2h. This property is well Illustrated In 
Figures S and 8 where the orbit separation of 1000 km gives the smallest 
baseline error. For the long baseline ranging, the baseline error is shown 
by Inequality (54) to increase directly with the orbit altitude and Inver- 
sely with the orbit separation. These properties are illustrated in 
Figures 5 through 13. In all cases, the range measurements must be taken 
at two widely separated satellite orbits in order to strengthen the 
geometry sufficiently to recover Che stations and the baseline at the cen- 
timeter level 1 10). 


The dependence of baseline error on the minimum elevation angle cannot 
be extracted directly from Figures 5 through 13 because the normalization 
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Figure 5. Nonulized baseline error versus baseline for different satellite 
orbit separations In the uncorrelated-path laser ranging systea. 
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Figure 6. Normalized baseline error versus baseline for different satellite 
orbit separations in the uncorrelated-path laser ranging system. 
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Figure 7. Noraaiized baseline error versus baseline for different satellite 
orbit separations In the uncorrelated-path laser ranging systea. 
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Figure 8. Nomalized baseline error versus baseline for different satellite 
orbit separations in the uncorrelated-path laser ranging system. 
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Figure 9. Nonielixed bescllne error vetsue besellne for different eacelilte 
orbit separation* in the uncorralated-path laser ranging systea# 
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Flgur« 10 * Norullz^d bMellnc error versus besellne for different setellite 
orbit seperstions in the uncorreleted**peth leser renglng systee* 
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Figure 11 


Noraallced bMelin* error vereuc beeeline for dlffereac eecellite 
orbic sepereclone in the uncorreleted-pech laeer ranging systea. 
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riff ir* 1? Normalized bazellne error verzue beeeline for different eetelllte 
* * orbit zeperetlone la the uocorreleted-peth laser ranging system. 
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Figure 13, NormeUsed baseline error versus baseline for different satellite 
orbit separations In the uncorrelated-path laser ranging system. 
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factor used In Equation (SI) Involvaa tha alnlMUi alavatlon angla which 

varlaa for dlffarant aata of paraaatara. In ordar to Invaatlgata this 

ralatlonahlp, tha abaoluta raa basallna arror la pLottad varaua tha 

baaallna dlatanca for varloua alnlaua alavatlon anglas In 

Figure 14; and varaua la plotted In Figure IS for varloua 

baaallna dlatancaa. Thaaa plota Indicate that the baaallna arror la a 

dacraaalng function of tha ainlauB alavatlon angla. Thla ralatlonahlp be*’ 

tween the baaellne error and the alnlaua elevation angle can be explained by 

the fact that while keeping the other paraaetera conatanc, the Increaae In 

alnlaua elevation angle will result In the shorter satelllte-to-ground- 

statlon distance. Consequently, the range aeasureaent errors, and hence 

the baseline error, will be saaller. 

In aost cases the orbit altitude and the alnlaua elevation angle are 

usually defined prior to the actual ranging, so the baseline error depends 

very auch on the orbit separation for a given baseline ranging. For 

exaaple, the altitudes for aany of the satellite orol 'i are 1000 ka. For a 

alnlaua elevation angle of 20* and an angle y of 60*, the lOOO-ka baseline 

has a raa error of around O.S <j - for a 1000-ka orbit separation and a raa 

AK 

error of around S for a 100-ka orbit separation (refer to Figure 2). 

This effect of orbit separation on the baseline error will even be aore 
pronounced for a longer baseline. 

4,3 Baseline Error for the Correlated-Path Laser Ranging 

If the weather conditions In both ground stations during the second 
satellite pass are alaost the saae as Chose during the first satellite 
pass, we would expect Che ataospherlc refraction Co be highly correlated. 
This Implies that 


BASELINE ERROR (<r, 
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Figure 14. Baseline error versus baseline for different elevation angles in 
the uncorrelated-path laser ranging system. 


BASELINE ERROR (or. 


50(J/a 


ORIGINAL PAGE IS 
OF POOR QUALITY 


46(Tfl 


420-O 


h«500km 
/* 100 Km 

x«o« 


38O3 


Bl *1500 km 


34^5 


30<To 


26o-fl 


Bl*IOOO km 


220fl 


I40g 


Bj_*500 km 


^5 10 15 20 25 30 35 40 45 50 55 60 

MINIMUM ELEVATION ANGLE. E„,„ (DEGREES) 


Figure 15. Besellne error versus minliaua elevaclon aogle for different baselines 
In the uncorrelated-*path Laser ranging systen. 
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p > 1 ffl ■ 1,2 . (55) 

Itt 

We refer to this cooflguraelon es Che correlaced-peCh ranging. The abso- 
lute and nonsalized ras baseline errors for this configuration can be 
obtained by subscltutlng Equation (55) into Equations (47) and (48). This 
yields, respectively. 


Og (corr) • 

Is 


^ <Jg cos Y 


■^V4(Bj^ cos Y + 1 )^)-Cj^ 2 l“ ■*‘V4(Bj^ cos Y * *)^0 

(56) 


and 


(corr) 


/2 cos Y sin B 


nln 


“AR 


Ih 


{Cii(h^ + ^ 


- + V4(Bj^ cos Y - 1)^1} 


(57) 


By substituting Equation (55) into Equation (49), we find Chat the minlaua 
baseline error of zero can be achieved when B, *0. 

By referring to Equation (53), the upper bound for the baseline error 
In this systea can be found as 
B, JT u ^ B, cos'^ Y sin E . 

^‘“">< 34 ^ 1 ^ ^ • »»> 

AR aln 


Since B^^^ increases with the orbit altitude (h) and decreases with the 
orbit separation (1), Inequality (58) Indicates chat the baseline error In 
this case also increases directly with the orbit altitude and inversely 
with Che orbit separation for long baselines. 

Equation (56) Is plotted versus Che baseline distance for four dif- 
ferent alniaua elevation angles in Figure 16; and it is plotted versus 


BASELINE ERROR (<r, 
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th« •inlaua •lavation «agl« In Flgura 17 to illustrnt* th« dapondanca of 
basallna arror on alnlaua alavation angla* Equation (57) la plottad varaua 
tha baaallna dlatanca In Flgurat 18 through 21 for dlffarant aats of 
orbit saparatlona, orbit altltudast ainlaua alavation anglaa and tha y 
anglaa. Thaaa plots Indlcata that tha ganaral propartlaa of tha basallna 
arror aantlonad In Sactlon 4.2 hava also hald trua hara for tha corralatad> 
path ranging. This Is to ba axpactad bacausa tha gaoaatrlc nodal usad In 
this sactlon, on which thasa propartlas ara basad. Is tha sana as that usad 
In Sactlon 4.2. 

It Is Intarastlng to nota tha laportanca of ataospharlc rafractlon 
error corralatlon In the satelllta laser ranging process. By referring to 
Equation (47), wa see chat Che higher Che correlation between Che refrac- 
tion errors on the two passes, Che lower Che baseline error. In par- 
ticular, Che baseline error for the uncorrelated-path ranging can be nore 
chan an order of nagnlcude higher chan chat for the correlaced-pach 
ranging. Since cha correlation batwaen rafractlon errors on cha two passes 
depends very auch on vO* ceaporal sapsraclon between chase passes. It Is 
laportanC Co keep Ch.^- '>ecwaan passes short enough In order to enhance 

Che refracclon arror c .elation. 


BASELINE ERROR 
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Figure 17* 


Baseline error versus minimum elevation angle for different baselines 
in the correlated’-path laser ranging system* 






Figure 18. Noraallred baseline error versus baseline for different satellite 
orbit separations in the correlated-path laser ranging system. 
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Figure 19. Noraelized baseline error versus baseline for different satellite 
orbit separations in the correlated-path laser ranging system. 
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S. COOROItlATE BRRORS FOR A SINGLE GROUND STATION 
In chls Mccion, w« •xaaloc tha coordioaea arrora for a aingla acacioo dua 
Co cha af facta of acaoapharic rafraccion and aaa how chay vary with raapact to 
tha orbital paraiaocara. By rafarring to cha coordinata covarianca aatrix of 
Equation (31), wa aaa that cha x-eoordinata varianca ia aaro, and tha y~ and 
z-coordinaca variancaa of a aingla ground atation ara givan by Equatlona (32) 
and (34), raspactivaly. Tha gaoaatry uaad to avaluata chaaa two variancaa ia 
alaoac cha aaiaa aa chat daacrlbad in Saction 4 (Ftgura 4), with cha axcaptiona 
that cha baaalina diatanca (B^^) batwaan two atatlona ia now raplaced by cha 
parpandlcular diatanca batriaan Cha x-axia and a ground atatlcn (Fq)* Thia 
gaoiaatry Is illuacratad in Figure 22. 

Tha foraulaa for cha varioua alavatlon anglaa given by Equatlona (41) to 
(45) are still valid hare with cha term (B^ cos y) replaced by 2/^. By 
aubacitutlng chase new equations into Equations (32), (34), (38) and (39), we 
can express cha coordinate variances aa 
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( 60 ) 


where 
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Plgur. 22. lUnging g«^««try for « .lngl« ground station. 
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U C. 


I •!« 2*^ + I 


■k‘1" 


‘4 


n - 1,?, 


(61) 


co« 


-I 


sin E 


sin E‘J 


n - 1.2, (62) 


n denotes the nth satellite pass, and p Is the g-error correlation coef- 
ficient associated with the measureaents acquired by the ground station 
during the two passes. 

By taking the square roots on both sides of Equations (59) and (60) and 
nonnalizing the results by the oaxlmua atmospheric refraction error 
we have the normalized rms coordinate errors 
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(63) 
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The approximate lower bounds for these coordinate errors can be found 


by assuming a zero y^ and letting be unity. This yields 
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(65) 


( 66 ) 


By coaparlag Equation (63) with Equation (48) » we can see that the charac*- 
terlatlcs of the y-coordlnate error follow exactly the characterlatlca of the 
baseline error. This Is to be expected because the y-coordlnate error Is the 
only error source contributing to the baseline error In the previous analysis. 
Equations (64) and (66) show that the z-coordinate error Is directly related 
to the satellite orbit separation (1) and is inversely related to the orbit 
altitude (h) for the small y^. Furthermore, Equation (65) indicates that a 
zero y-coordinate error can only be obtained when the refraction errors during 
the two passes are positively correlated (p • +1); whereas. Equation (66) 
indicates that the negatively correlated refraction errors (p • -1) are 
ne». assary for the zero z-coordinate error. These results can be explained by 
the geometry of the ranging system. When y^^ Is zero, the satellite ground 
tracks are equidistant from each side of the ground station. If the refrac- 
tion errors are positively correlated, this will give rise to the y-coordlnate 
range errors of equal magnitudes and opposite directions on the two passes, 
and therefore, the resultant y-coordlnate error for the station will be zero. 
Since both satellite orbits are above the station at equal altitudes, in order 
to have a zero z-coordinate error for the station (which means the z- 
coordinate range errors on the two passes are equal in magnitude but opposite 
in direction), the refraction errors during the two passes must be negatively 


correlated. 
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When jTq Is large compared to t, the coordinate errors can be approximated 
as 
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3 sin E 


min 
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T 

\ 



(cj - ZpCjCj + C^) 




(67) 


and 


3 sin E 


min 


'AR 



(cj - ZpCjC^ + c\) 


Itl 


( 68 ) 


Equation (67) indicates that the y-coordinate error Is proportional to the 
square of y^. Whereas, Equation (68) indicates that the z-coordinate error is 
a linear and cubic functions of y^. By comparing these two equations, we find 
chat the z-coordinate error for a ground station can be a lot larger chan the 
y-coordinate error if Che satellite altitude (h) is small compared to y^^. 

When h and y^ are comparable, Che coordinate errors of equal magnitudes are 
resulted. In this case, both coordinate errors are related directly to the 
satellite altitude (h) and inversely to the orbit separation (1). It should 
also be noted Chat when y^ is large, and are approximately equal, and 
Che minimum coordinate errors can be obtained when p equals unity. This 
means Chat the smallest possible station coordinate errors for large y^ can be 
achieved when the atmospheric refraction errors are positively correlated 
( this is analogous to the result obtained in Section <* where the correlated- 
path ranging always gives the smallest baseline error). 
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6. CONCLUSION 

BM«d upon Ch« standard regrasslon oodel and tha gaoaatry of the 
ranging sites, the general expression for the baseline error In the spheri- 
cally syonetrlc ataoaphere la derived. 

It has been shown that the y angle and, hence, the angle of Intersection 
between the baseline and the satellite ground tracks play an Important 
role In the baseline error determination. In particular, the minimum 
baseline error can be achieved when y approaches 90*, that Is, when the 
baseline or the two ground stations Is parallel to the satellite ground 
tracks. Besides the fact that the baseline error has a strong dependence 
on the angle y, the choices of minimum elevation angle, orbit altitude and 
orbit separation are also Important In determining the magnitude of this 
error. 

Figures (5) through (26) indicate that the baseline error grows rapidly 
as the baseline increases, and it can be more than an order of magnitude 
larger than the refraction error for the intercontinental baselines. In 
cases where the baselines are short, the baseline errors can be a fraction 
of the atmospheric refraction errors. 

The analysis In this report does not consider the baseline error for 
the nonhomogeneous atmosphere. In the nonhomogeneous atmospheres, addi- 
tional refraction error terms due to the gradient correction must be 
Included. These terms may have a significant contribution to the total 


baseline error. 
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